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Abstract: The transverse momentum transfer correlation is introduced as a sensi- 
tive probe that can be used to discrimate between models for parton dynamics in 
low-x deep-inelastic scattering. Expectations for uncorrelated models and models 
ON . with short-range or long-range correlations are discussed and confronted to results 

^ ! obtained from the Lepto and Ariadne Monte Carlo simulation programmes. 
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1 Introduction 

^ : 

The successful description of the nucleon structure function by perturbative QCD, using the 
DGLAP (Dokshitzer-Gribov-Lipatov-Altarelli-Parisi) parton evolution equations constitutes 
^ . a major success of QCD. However, at very small Bjorken-x, these equations are expected to 
>- ! become invalid. An alternative ansatz for the small-a: regime is the BFKL (Balitsky-Fadin- 
Kuraev-Lipatov) equation 0. At lowest order, the BFKL and DGLAP equations resum the 
leading logarithmic (a s lnl/x) n ) or (a s \nQ 2 /Q 2 ] ) n ) contributions, respectively, with Q 2 being 
the virtuality of the exchanged photon in deep-inelastic neutral current ep collisions. The 
leading-log DGLAP ansatz implies strong ordering (Ql <C kj\ ■ • • <C kxl <C - • • Q 2 ) of the 
transverse momenta (k^i in the parton cascade. Here, (/ct)i is the transverse momentum 
(measured relative to the incident proton direction) of the i-th emitted "final state" parton (see 
Fig. [I]). According to BFKL kinematics, the transverse momenta are no longer strongly ordered 
but may be thought of as forming, roughly speaking, a random walk with (fcf^j ~ (kj)i+\. 

Measurements of the hadronic final state resulting from the initial-state parton cascade 
should be sensitive to the type of evolution. For example, without strong fc^-ordering, more 
transverse energy is to be expected from the BFKL than from the DGLAP type of evolution 
in a rapidity interval between the proton fragmentation region and the current fragmentation 
region. A similar expectation holds for the inclusive transverse momentum distribution of 
hadrons (0,3- So far, however, none of the studied observables (including recent measurements 
of forward jet and large-pr charged and neutral pions ||) are directly probing the correlation 
structure of the parton cascade. To allow additional discrimination between a strictly kx- 
ordered and an unordered scenario, it is necessary to go beyond single-particle observables and 
consider directly measures of the transverse momentum correlations in the parton cascade. 
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Figure 1: Parton evolution in the ladder approximation. The longitudinal fractional momenta 
Xi and transverse momenta kn of subsequently emitted partons are indicated. 

2 The transverse momentum transfer 

In this paper we consider the transverse momentum transfer across rapidity y, H(y)- This 
quantity is defined as the vector-sum of the transverse momenta fcj of all particles with rapidities 
smaller than y in an event with n particles. For simplicity, we consider in the following only 
one of the components of the transverse momentum vectors (denoted by hi) and define 

n 

r%) = X>%-^)> (!) 

i=l 

where 8(x) = for x < and 8(x) = 1 for x > 0. H(y) is a random function, varying from 
event to event. A set of k-th order moment-functions can be constructed defined asQ 

(n( yi )n(y 2 )...n(y fe )>, (2 ) 

where () means that the average is taken over an ensemble of events. From the moment- 
functions, one can further define a set of correlation functions by taking a cluster decomposition 
(see eg. 0). The two-point correlation function D 2 (y 1 ,y 2 ) is then given by 

D 2 (y u y 2 ) = (n( yi )n{y 2 )) - (U( yi )) (I% 2 )). (3) 

In particular, for yi = y 2 = y we define D 2 {y) = D 2 (y,y). This quantity measures the 
variance of the transverse momentum transfer distribution across the rapidity "boundary" y. 
D 2 (y) provides a measure of the rapidity structure of the correlations between the transverse 
momenta as can be seen from the relation || 

/y /"Go f 

dyi dy 2 d 2 k 1 d 2 k 2 k 1 ■ k 2 P2(h,yi, k 2 ,y 2 ), (4) 
-oo J y J 

where p 2 (ki,yi, k 2 ,y 2 ) is the inclusive two-particle density. 

With reference to the parton- level diagram of Fig. [I], it should be noted that H(y) is in that 
case equal to the (propagator) transverse momentum exchanged between the vertices % and 
i + 1 if the position y is located between particles i and % + 1 on the rapidity axis. Transverse 

^ee Q for a discussion of the theory of random functions 
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momentum transfer correlations are, therefore, a direct measure of the correlations between the 
exchanges. 

In models of uncorrelated production, the transverse momentum of a given particle can be 
compensated by that of any other particle or group of particles. Consequently, the correlation 
length describing the transverse momentum correlations, which follow from momentum conser- 
vation alone, is of the order of half the available rapidity range and should thus increase with 
energy or W in DIS. For uncorrelated transverse momenta D 2 (y) is expected to increase when 
y moves from the edge of the rapidity space to y — since, for independent random variables, 
the variance of a sum increases proportionally to the number of terms included. For the same 
reason also D 2 (0) should increase with the number of particles produced. 

Going beyond the uncorrelated production assumption, one could assume that the relevant 
part of the diagram in Fig. [I] (neglecting any dependence on longitudinal momenta) can be 
schematically written as a product of nearest-neighbour correlated "links" 

|x„| 2 ^n v(u i ,n i+1 ), (5) 

i 

where the functions V represent the (iterated) kernel of the diagram. In that case, an important 
role is played by the eigenvalues and eigenf unctions of the kernel, considered as a transfer 
integral operator ||. For a sufficiently long ladder, it can be shown that the function D 2 (yi, y 2 ) 
is controlled by the largest eigenvalue Aoi leading to a behaviour 

D 2 ( Vl , y 2 ) « exp {-p(l - A i)||/i -y 2 \}, (6) 

with p the mean parton multiplicity per unit of rapidity. In this short-range model, the corre- 
lation length is therefore given by l/p(l — Aoi). 

In models with short-range correlations, by definition, particles far away in rapidity cannot 
be correlated. As a result, the compensation of transverse momenta is local, i.e. it is compen- 
sated by its neighbours. The corresponding correlation length A is energy independent. These 
general properties are reflected in the behaviour of D 2 (y). Only particles in the rapidity interval 
A = (—A + y,y + A) contribute to the transverse momentum transfer at y. Therefore, D 2 (y) 
is expected not to increase faster with energy than the particle density in the interval A (for 
constant A) which is known to increase only logarithmically. 

In a ^-ordered cascade, correlations will be long-range. The correlation length can be 
expected to increase with the 'length of the cascade' in rapidity space, i.e. proportional to 
lnl/x or ha.W (at fixed Q 2 ) in DIS. For the case of the BFKL scenario, it can further be 
expected that the correlation length is related to the dominant eigenvalues of the BFKL kernel 
(cfr. Eq. |) 

3 A Monte Carlo study 

Predictions for the cases of the ordered and unordered cascades are obtained from DIS Monte 
Carlo models, which incorporate the QCD evolution in different approximations and utilize 
phenomenological models for the non-perturbative hadronization phase. The MEPS model 
(Matrix Element plus Parton Shower) [IIJ, incorporates the QCD matrix elements up to first 
order, with additional soft emissions generated by adding leading log parton showers. In the 
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colour dipole model (CDM) JTT , 12 ] radiation stems from colour dipoles formed by the colour 



charges. Both programs use the Lund string model |13| for hadronization. The CDM description 



of gluon emission is similar to that of the BFKL evolution, because the gluons emitted by the 



dipoles do not obey strong kx ordering fl4 |. In MEPS the partons are strongly ordered in kx, 
because they are based upon leading log DGLAP parton showers. 

Two samples of one million events have been generated using the latest versions of the 
models (Lepto 6.5 for MEPS, Ariadne 4.08 for CDM) with the parton density parametrisation 
GRV-94 [[HJ. The fraction of events undergoing a Soft Colour Interaction |16[] was set to zero. In 
all figures the transverse momentum transfer was calculated using final state hadrons, including 
neutrals so that transverse momentum is conserved for the full event. 

Figure ^ shows the variance of the transverse momentum transfer, D 2 (y), scaled to the 
average hadron transverse momentum squared in different bins of Q 2 and x. Large differences 
between models are observed at low x. Because the variance is scaled, this effect is not due 
to the larger p^-flow produced by Ariadne. As expected, the transverse momentum transfer 
fluctuations are larger in the unordered scenario. 

The correlation function, computed as the ratio D 2 (yi,y2)/D(yi)D(y 2 ), is shown in Fig. |^. 
The correlation length (shown in Fig. |j) is calculated along slices perpendicular to the diagonal 
Hi = U2- Within these slices the correlation decreases exponentially with \yi — y 2 \ and the 
correlation length is defined as the slope of an expontial fit. The models used in this simulation 
show no strong dependence of the correlation length on x or Q 2 . Again, differences are observed 
in the behaviour of the correlation length between the Lepto and Ariadne Monte Carlo models. 
In particular, Ariadne predicts a symmetric behaviour in the proton and photon fragmentation 
hemisphere with a double-peaked structure for the correlation length, while Lepto predicts 
longer correlation lengths in the photon fragmentation hemisphere. 



4 Conclusion 

It has been demonstrated that the transverse momentum transfer correlation is a theoretically 
attractive variable to discriminate between models for QCD evolution in parton cascades. The 
two considered models, Lepton and Ariadne, show large differences in the behaviour of the 
transverse momentum transfer fluctuation and correlation length. However, further work is 
needed to investigate the feasibility of an experimental measurement. 
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Figure 3: The fluctuation of the transverse momentum transfer, D 2 , scaled to the average 
transverse momentum squared, (pj<), is shown as a function of rapidity y in different bins of 
Q 2 and x for Lepto (full line) and Ariadne (dashed line). 
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Figure 4: The correlation length, A, defined as in the text, is shown as a function of rapidity y 
in different bins of Q 2 and x for Lepto (full line) and Ariadne (dashed line). 



7 



